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A SIMPLE FINITE DIFFERENCE APPROACH USING.UNSTRUCTURED MESHES 
FROM FEM MESH GENERATORS 


F. Anibal Fernandez*, Lukasz Kulast 


Abstract: This paper presents a methodology for the use of Finite Differences in two-dimensional eigenproblems 
using a mesh, of triangles produced for the finite element method. Expressions are given for the first and second 
derivatives at a grid point in terms of nedal values at neighboring points. The method has been tested analysing a 
waveguide resonator including a dielectric rod. 


1. Introduction 


In all methods involving domain discretization, the accuracy and quality of the solution depend crucially 
on the mesh size. However, in most problems the degree of discretization required to achieve a certain 
level of accuracy varies across the problem domain. Therefore, the efficiency of the method relics on its 
ability of using meshes that can be adapted to the problem. The Finite Element Method (FEM) has 
become one of the most powerful numerical tools in engineering and physics because of its flexibility and 
versatility, and this is partly due to its ability of using adaptive meshes, but has a serious limitation in 
efficiency for large number of variables. Standard Finite Difference methods offer a simpler formulation 
and implementation than finite element methods, but they require a uniform mesh. This is inappropriate 
in cases where the solution sought varies widely in the problem domain or for structures containing small 
and/or odd-shaped features or with boundaries non-conforming to the grid. Subgridding [2], or local grid 
refinement, [1] has been used in finite difference methods. Nonuniform and unstructured meshes [1] are 
also possible and they provide the highest flexibility in adapting mesh to a given problem, which has 
proven its usefulness in many practical applications [{1, 3]. 

Among unstructured meshes used in FD triangular are often preferable, as they can be generated 
using existing FEM mesh gencrators. Existing schemes using triangular or tetrahedron meshes are often 
complicated as they iterate unstructured and uniform meshes separately [3]. In this work we have used a 
meshing method developed for Finite Element calculations in 2D, which allows the construction of meshes 
adapted to a given density function or if included in a calculation process, to produce an adapted mesh 
(4]. Expressions for the first and second order partial derivatives of a function at a mesh node have been 
derived in terms of the function values at a number of neighboring mesh points for the implementation of a 
finite difference method. The method reduce to a regular FD scheme, when applied to a uniform grid, thus 
no spatial interpolation has to be performed between unstructured and uniform meshes and expressions 
for implementing Neumann boundary conditions, periodicity conditions and absorbing conditions can 
also be added as they exist for a regular FD scheme [1]. The procedure have been tested by calculating 
resonant frequencies of an empty cavity and one containing a dielectric rod, comparing these results with 
those from a traditional finite differences approach. 


2. The Theory 


Considering. first briefly the one-dimensional case, we have that the expressions for the first and second 
derivatives that arise from fitting a second order curve to three consecutive points are: 


fis { [22;(a; — G41) + (02,1) + (03 - 3) + fina 
+ [2ei (aig, — @i1) + (a2, — 22,)] - fi 
+ [2aj(ay-1 ~ 2i) + (2? — 2?_,)] + fig Mia (1) 
Al = 2 (wi — wins) + fina + (Wigs — Bin) - fi + (Win — 24) - fit) /A 
(2) 
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which reduce to the central difference expressions for the first and second derivatives. The same idea can 
be applied to the two-dimensional case. Here we start with a mesh of triangles designed for finite element 
calculations and fit a second order surface to each mesh node. Since a set of six points is needed in order 
to fit a unique surface, 5 immediate neighboring nodes can be selected using the connectivity information 
provided by the mesh generator. Defining the surface f by: 


f(x,y) = ag2? + asy” + aary + 3% + ary + a4 (3) 
the resultant equations can be written as: 
f = Ma (4) 


where the vectors f and a are the function values at the six nodes, and the coefficients in (3), respectively. 
The matrix M of order six can be inverted analytically finding the values of the required coefficients. 
With this, expressions for the first and second order derivatives can be found: 
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These expressions also reduce to the standard formulae when applied to a regular mesh. This allows one to 
easily adapt the presented algorithm to a case, when only part of the domain is covered by an unstructured 
triangular mesh. Although this extension is not considered here, its application is straightforward and 
does not require spatial interpolations at the interface between unstructured and regular meshes. 


3. Implementation in 2D FDFD eigenproblem formulation 


To implement the presented algorithm into the 2D FD-eigenproblem formulation let us consider a 2D 
rectangular resonator, as presented in Fig. 1. If we assume a TMjnrno Mode, we can formulate a following 
equation describing FE, field component along x and y directions 


oF OF 


2 2 _ 
VEL —k B.= (sat 5a 


) E, —w*yeE, =0 (6) 
After the mesh is generated, which for the considered structure is presented in Fig. 3, for each mesh point 
one should find 5 neighboring points and solve numerically or analytically system (4) to obtain a@)-a¢. 
Using expressions describing derivatives (5) equation (6) written for each point form the following set of 
equations: 

Ae— w* pgéoDe =0 (7) 


where matrix A and diagonal matrix D, contain coefficients a5, ag and the effective permittivity cal- 
culated for each mesh point. A local value of é¢fs may be easily calculated using data provided by a 
mesh generator - a value of permittivity inside each triangle. System (7) containing unstructured mesh 
operator A can now be solved using a standard eigenvalue solver. 


4. Results 


Two examples are presented here to demonstrate the validity of this approach. First, an empty cavity 
is analysed both using a standard FD discretization and the method proposed here. For TM,»9 modes 
the problem is solved in 2D. The cross-section dimensions of the cavity are 6mm x 5mm. The solutions 
obtained are presented in Table 1. This problem is particularly well suited for the standard FD method 


Regular grid 


587 


Triangular grid 


Number of grid points | error [%] 
0.36 
0.04 
0.32 
1167 0.03 


Table 1: Errors in computations of the resonant frequency for the TMiio mode obtained using Finite Differences for regular and 


unstructured meshes. 


and despite this, the proposed method shows more accurate results than in case of a regular mesh having 


13 x 11 grid points. 


The second example consists of the same cavity but with a full-height dielectric rod of lmm x Imm 
square cross-section and relative permittivity 8 inserted at the center of the cavity and rotated 45 degrees 
with respect to the walls. Fig.1 shows a diagram of the structure. The field distribution for the FE, 


Fig. 1: Geometry of the loaded cavity 


component of the lowest TM mode is shown in Fig.2 and Fig.3 shows the non-uniform grid used in our 


calculations. For comparison, results were 
108661 grid points (361 x 301 grid nodes). 


calculated with a uniform grid of increasing size, of up to 
The resonant frequency from the unstructured grid, which 


uses 301] points is 26.401G Hz, corresponds to an error of 0.0335% compared to the results obtained with 
the finest uniform grid. This error is equivalent to that of a uniform grid of more than 15000 points 


(145 x 121 grid nodes). 
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Fig. 2: Ez, component of the lowest TM mode. 
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Fig. 3: 2D grid used in the calculation containing 301 nodes 
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5. Conclusions 


A 2D algorithm has been presented for the implementation of the finite difference method using unstruc- 
tured grids designed originally for finite element calculations, which can be easily adapted to complicated 
problem geometries. Results show considerable advantages over the standard FD method, when applied 
to a problem not conforming to a regular grid. The methodology can easily be adopted to a regular FD 
algorithm having unstructured mesh only locally. 
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